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I. INTRODUCTION 



The development of techniques for cooling and trap- 
ping of a wide variety of atomic and molecular species 
in recent years has created exciting opportunities for 
probing and controlling atomic and molecular encoun- 
ters with unprecedented precision While many of 
the initial studies of cold atoms and molecules were 
centered on the creation of dense samples of cold and 
ultracold matter, more recent work has focused on the 
manipulation and control of intermolecular interactions. 



with the ultimate aim of achieving quantum control of 
atomic and molecular collisions [2]. Though the ideas 
of quantum control of chemical reactions were proposed 
many years ago, the ability to create ultracold molecules 
in specific quantum states has given further stimulus 
to this field. Its development requires that molecular 
properties and collisional behavior be well understood 
at cold and ultracold temperatures where the dynamics 
of molecules are dramatically different compared to col- 
lisions at elevated temperatures. Over the last ten years 
significant progress has been achieved both in theoretical 
and experimental works. The experimental methods 
such as photoassociation spectroscopy, magnetic tuning 
of Feshbach resonances, buffer-gas cooling, and Stark 
deceleration have been developed and applied to a 
variety of molecular systems. Novel methods to study 
ultracold chemical reactions involving ion - molecule 
systems in a linear Paul trap have been proposed Q. 
External control of chemical reactions using electric and 
magnetic fields is another area of active interest Q. 
The aim of this chapter is to provide an overview of 
recent progress in characterizing molecular processes and 
chemical reactions at cold and ultracold temperatures 
with particular emphasis on theoretical developments 
in quantum dynamics simulations of atom - molecule 
collision systems over the last ten years. 

In contrast to scattering at thermal energies, ultracold 
collisions offer fascinating and unique opportunities to 
study molecular encounters in the extreme quantum 
regime where the entire collision can be dominated by 
a single partial wave. One of the main motivations of 
current experimental efforts to create dense samples of 
ultracold molecules is to study the possibility of chemical 
reactions at temperatures close to absolute zero. While 
Wigner's law [3, [lOj predicts that rate coefficients of 
exothermic processes are finite in the zero-energy limit 
it does not say if the rate coefficient will be large enough 
for reactions to be observable in an experiment. Nor does 
it say anything about how the rate coefficients at zero 
temperature depend on the interaction potential. Most 
chemical reactions between neutral atoms and molecules 
involve an energy barrier and it is not clear if chemical 
reactions between them generally occur with measurable 
rates at ultracold temperatures. Calculations for the 
F -|- H2 reaction which proceeds by tunneling at low 
energies have shown that the reaction may occur with a 
significant rate at ultracold temperatures [Tl|. There is 
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also experimental and theoretical evidence that chemical 
reactions involving heavy-atom tunneling of carbon [T2j 
and fluorine atoms can occur with significant rate 
coefficients at low temperatures. 

Photoassociation experiments involving alkali-metal 
systems have stimulated considerable interest in chem- 
ical reactivity in alkali-metal dimer - alkali-metal atom 
collisions at ultracold temperatures . Rearrangement 
collisions in identical particle alkali-metal trimer systems 
occur without energy barrier and recent studies have 
demonstrated that chemical reactions in alkali-metal 
atom - alkali-metal dimer collisions may be very fast 
at ultracold temperatures. Unlike tunneling-dominated 
reactions, the limiting values of the rate coefficients 
for alkali-metal systems are less sensitive to vibrational 
excitation of the dimer. 

In this chapter we give an overview of recent studies 
of ultracold atom - molecule collisions focusing on non- 
reactive and reactive systems and the effect of vibrational 
excitation of the molecule on the collisional outcome. We 
will discuss both tunneling-dominated and barrierless re- 
actions and examine recent efforts in extending these 
studies to ionic systems as well as molecule - molecule 
systems. We consider mostly the novel aspects of colli- 
sional dynamics of atom - diatom systems at cold and ul- 
tracold temperatures with illustrative results for specific 
systems. For more comprehensive discussion of cold and 
ultracold collisions including reactive and non-reactive 
processes and the effect of external fields we refer the 
reader to several review articles [1, [l3l - [l5| that have 
appeared in the last few years. For details of the theoret- 
ical formalisms we refer to the chapters by Hutson and 
by Tscherbul and Krems. 



II. INELASTIC ATOM - MOLECULE 
COLLISIONS 

Theoretical studies of ultracold molecules received in- 
tense interest after the success of ultracold atom pho- 
toassociation and buffer gas cooling experiments which 
demonstrated that a wide array of molecular systems 
with thermal and non-thermal vibrational energy dis- 
tributions can be created in ultracold traps. The col- 
lisional loss of trapped molecules is an important issue in 
these experiments. Photoassociation produces molecules 
in highly excited vibrational levels. Whether the excited 
molecules decay by vibrational quenching or through 
chemical reaction is an intriguing question. While an 
extensive literature exists on the collisional relaxation 
of vibrationally excited molecules at elevated tempera- 
tures not much was known on the magnitude of the re- 
laxation rate coefficients at temperatures lower than one 
Kelvin. Though a few earlier reports jT6l - [l8| on atom 
- diatom collisions in the Wigner threshold regime had 
been published, a detailed investigation of the depen- 



dence of the relaxation rate coefficients on the internal 
energy of molecules and their sensitivity to details of the 
interaction potential has not been carried out. Here, we 
give a brief account of recent quantum dynamics calcu- 
lations of vibrational and rotational energy transfer in 
atom - diatom collisions at cold and ultracold tempera- 
tures. We focus on a few representative systems to illus- 
trate the main features of energy transfer in non-reactive 
atom - molecule collisions at ultracold temperatures and 
we show how the corresponding rate coefficients are in- 
fluenced by rotational or vibrational excitation of the 
molecule. 



A. Vibrational and rotational relaxation 

1. Collisions at cold and ultracold temperatures 

As discussed in the chapter by Hutson, at very low 
energies, scattering is dominated by s- waves and the 
scattering cross section can be expressed in terms of 
a single parameter called the scattering length. For 
single-channel scattering where only elastic scattering 
is possible, the scattering length is a real quantity 
and the magnitude of the cross section in the s-wave 
limit is given by cr = Aira'^ where a is the scattering 
length. For multichannel scattering, as in vibrationally 
or rotationally inelastic collisions of molecules, the 
scattering length is a complex number and it is denoted 
as a^j = a-uj — iPvj where v and j are, respectively, the 
initial vibrational and rotational quantum numbers of 
the molecule [l^, [l^- The limiting value of the elastic 
cross section in the presence of inelastic scattering is 
given by al^^ = 47r|a„jf = 47r(a2^. + Pl^). The total 
inelastic quenching cross section from a given initial 
rovibrational level of the molecule is related to the 
imaginary part of the scattering length through the 
relation cr*" = 47r/3i,j/k„j where k^j is the wave vector 
in the incident channel. The quenching rate coefficient 
becomes constant at ultralow temperatures and is given 
by fc^j = 47r/i/3t,j//j, at zero temperature where /i is 
the reduced mass of the collisional system. Thus, the 
rovibrational relaxation rate coefficients attain finite 
values for different initial vibrational and rotational 
levels of the molecule. The dependence of the rate 
coefficients on v and j has been an important issue in 
cold molecule research because exothermic vibrational 
and rotational relaxation collisions are a major pathway 
for trap loss in cooling and trapping experiments. 

Initial studies of rotational and vibrational relaxation 
of atom - molecule systems at cold and ultracold tem- 
peratures have mostly focused on van der Waals systems 
suchasHe-H2 [H-il, He-CO [Hill, andHe-Os 
Owing to the importance of some of these systems in as- 
trophysical environments, extensive calculations of low- 
temperature behavior of rate coefficients have been per- 
formed for coUisions of H2 [l9| and CO [H H with 
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both ^He and "^He. For both systems reasonably accu- 
rate intermolecular potentials have been reported. The 
initial calculations on the He— H2 system employed the 
potential energy surface (PES) of Muchnick and Russek 
(MR) For He-H2 vibrational excitation of the H2 

molecule has a dramatic effect on the zero-temperature 
quenching rate coefficients. As illustrated in Fig. [U the 
vibrational quenching rate coefficients increase by about 
three orders of magnitude between v = 1 and w = f of 
the H2 molecule Il9|. 




Temperature (K) 



FIG. 1: Rate coefficients for the quenching of H2(«, j = 0) by 
collisions with He atoms as functions of the temperature for 
1; = 1 — 10 of the H2 molecule. Reproduced with permission 
from Balakrishnan et al. [l9| . 



The quenching rate coefficients exhibit a minimum at 
around 10 K which roughly corresponds to the depth of 
the van der Waals interaction potential. This behavior 
appears to be a characteristic of vibrational quenching 
rate coefficients. For incident energies lower than the well 
depth the rate coefficient exhibits a minimum and with 
subsequent decrease in temperature the rate coefficient 
begins to increase before attaining the Wigner limit. For 
systems with deeper van der Waals wells the minimum is 
shifted to higher temperatures. Measurements of vibra- 
tional relaxation rate coefficients for the H2— CO system 
have confirmed this behavior 1271. 



Balakrishnan, Forrey, and Dalgarno [28| also investi- 
gated vibrational relaxation of H2 in collisions with H 
atoms for vibrational quantum numbers w = 1 — 12 of 
the H2 molecule. They adopted a non-reactive scattering 
formalism and neglected the rotational motion of the 
H2 molecule. The calculations showed that vibrational 
relaxation rate coefficients are strongly dependent on 
the initial vibrational level of the H2 molecule. The 
relaxation rate coefficients were found to increase by 
about seven orders of magnitude between vibrational 
levels V = 1 and v = 12. The dramatic variation in the 
rate coefficients with increase in vibrational excitation 
was explained in terms of the matrix elements of the 



interaction potential between the vibrational wavefunc- 
tions as functions of the atom - molecule center-of-mass 
separation. 

Vibrational relaxation rate coefficients in atom - 
molecule systems are often influenced by van der Waals 
complexes formed during the collision process. Decay 
of these complexes leads to resonances in the energy 
dependence of the relaxation cross sections (see Fig. [2] 
for Ar— D2 collisions). Applying effective range theory to 
describe ultracold collisions of the He— H2 system, Bal- 
akrishnan et al. [1^ and Forrey et al. [23| demonstrated 
that vibrational pre-dissociation lifetimes of resonances 
that lie close to the energy threshold can be derived 
accurately from the value of the zero-temperature 
quenching rate coefficient. This formalism was extended 
to describe vibrational relaxation of trapped molecules 
and it was shown that the vibrational relaxation rate is 
controlled by the most weakly bound state of the van 
der Waals complex (2^. In a related work Dashevskaya 
et al. 29] have shown that vibrational quenching of 
H2(w — 1, j = 0) at low temperatures can be described 
using a two channel approximation within the quasi- 
classical method provided appropriate parameters are 
employed in the calculations. Subsequently, Cote et 
al. [so] generalized this method to predict vibrational 
relaxation lifetimes of atom - diatom van der Waals 
complexes with energies near the dissociation threshold. 




Kinetic energy (cm ) 



FIG. 2: Cross sections for quenching of the « = 1, j = level 
of D2 in collisions with Ar atoms resolved into the different 
rotational levels j' in v' — as functions of the incident kinetic 
energy. Reproduced with permission from Uudus et al. ■ 



Unlike in thermal energy collisions, the presence 
of a weakly bound state in the vicinity of a channel 
threshold can dramatically influence the cross sections 
in the ultracold regime. This is illustrated in Fig. [5] 
where the cross sections for vibrational relaxation of 
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D2(u = l,j = 0) in collisions with Ar atoms jsTj are 
presented over an energy range of 10~^ — 10'^ cm~^. 
The cross sections exhibit a curvature characteristic of a 
resonant enhancement in the energy range 10~^ — 10"'^ 
cm~^. Such enhancement of the cross section can occur 
when the interaction potential supports a virtual state 
or a very weakly bound state near the channel threshold 
leading to a zero-energy resonance. The virtual state 
is characterized by a large negative scattering length 
while the bound state is characterized by a large positive 
scattering length. For the present case the real part of 
the scattering length is large and positive (aio = 97.0 
A) and the resonance occurs from the decay of a loosely 
bound van der Waals complex supported by the entrance 
channel potential. For energies below 10~^ cm~^ the 
cross section is dominated by s-wave scattering in the 
incident channel and the zero-energy resonance arises 
from s-wave scattering in the entrance channel. 

The resonant enhancement is more clearly seen in the 
plot of the reaction probability as a function of the kinetic 
energy shown in Fig. [3l The probability peaks at an en- 
ergy of 6.0 X lO"** cm^^ which roughly corresponds to the 
binding energy of the quasibound state. The resonance 
appears in the scattering calculations at energies above 
the threshold due to its close proximity to the channel 
threshold. The binding energy of the quasibound state 
can be estimated using the scattering length approxima- 
tion [l^, [3l|. The magnitude of the binding energy is 
given by \Eb\ = h? cos27io/(2/i|aio|^) where fi is the re- 
duced mass of the Ar— D2 system, aio ~ aio ~ i/3io is 
the scattering length for the v = l,j = level, and 
710 = tan~^ (/3io/q^io)- This yields a value of \Ei,\ — 
4.9 X 10"'' cm~^, in reasonable agreement with the exact 
value derived from scattering calculations. 

A more accurate value of the binding energy can be ob- 
tained using the effective range formula given by Forrey 
et al. [13]: 

where tq is the effective range of the potential which 
may be evaluated by fitting the low energy behavior of 
the phase shift for the elastic channel to the standard 
effective range formula, fcio cotJio — — 1/aio + 
For Ar— D2(w = l,j = 0) collisions, the effective range 
formula yields rg = 16.32 A. The resonance position 
calculated using the effective range approximation is 
\Ei,\ = 5.95 X 10~^ cm~^, in excellent agreement with the 
value of 6.0 x 10~^ cm~^ obtained from the scattering 
calculations. It is generally very difficult to accurately 
evaluate energies of such weakly bound states using 
standard bound state codes and the effective range 
formula provides a convenient and reliable method to 
calculate binding energies of weakly bound states that 
lead to zero-energy resonances. 
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FIG. 3: Total probability of quenching of the v = 1, j — level 
of D2 in collisions with Ar atoms as a function of the incident 
kinetic energy. The peak value of the probability corresponds 
to a zero energy resonance. Reproduced with permission from 
Uudus et al. [3ll |. 



One of the challenging aspects of cold and ultracold 
collisions is the sensitivity to details of the interaction 
potential. Even the best available methods for the 
electronic structure calculations of PESs result in errors 
much larger than the collision energies in the cold and 
ultracold regime and the dynamics calculations are often 
sensitive to small changes in the interaction potential. 
To explore the sensitivity of cold and ultracold collisions 
to details of the interaction potential, Lee et al. [s^l 
performed a comparative study of the ultracold collision 
dynamics of the He— II2 system using the MR potential 
and a more recent ab initio potential developed by 
Boothroyd, Martin and Peterson (BMP) [Hi. The 
BMP potential was considered to be an improvement, 
approaching chemical accuracy, over all conformations 
compared to the MR potential. However, significant 
differences were observed for vibrational relaxation of 
the V = l,j = state of the H2 molecule in collisions 
with He computed using the two surfaces. The limiting 
value of the quenching rate coefficient on the BMP 
surface was found to be about three orders of magnitude 
larger than that of the MR surface. The difference 
was attributed to the more anisotropic nature of the 
BMP surface leading to larger values of the off-diagonal 
elements responsible for driving vibrational transitions. 
Indeed, it was found that the vibrational quenching of 
the V = I, j = level was dominated by the transition to 
the v' = 0,j' = 8 level which is driven by the high-order 
anisotropic terms of the interaction potential. 

To explore the behavior of inelastic collisions involving 
polar molecules at ultracold temperatures, Balakrishnan, 
Forrey and Dalgarno ^23^ investigated vibrational and ro- 
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tational relaxation of CO in ^He— CO collisions. The dy- 
namics of the He— CO system was found to exhibit signif- 
icantly different features at low temperatures compared 
to the He— H2 system. Quantum scattering calculations 
of *He and "^He collisions with the CO molecule revealed 
that the larger reduced mass and the deeper van der 
Waals interaction potential of the He— CO system give 
rise to a number of shape resonances in the energy de- 
pendence of vibrational relaxation cross sections [23|, [13] . 
The effect of shape resonances on low temperature vi- 
brational relaxation rate coefficients will be discussed in 
the next subsection. The computed values of vibrational 
relaxation rate coefficients for both "^He and ''He colli- 
sions with CO(t; — 1) have been found to be in jpod 
agreement with experimental data of Reid et al. 34] in 
the temperature range 35 — 100 K. Calculations of vibra- 
tional relaxation rate coefficients for the He— CO system 
in the temperature range 35 — 1500 K have also been re- 
ported by Krems [35]- He has shown that inclusion of 
the centrifugal distortion of the vibrational wavefunction 
enhances the relaxation process, and that the quenching 
rate coefficients are sensitive to high-order anisotropic 
terms in the angular expansion of the interaction poten- 
tial [3^. 

Bodo, Gianturco and Dalgarno [13] have extended the 
work of Balakrishnan et al. [2^ to study the vibrational 
relaxation of excited C0(?; = 2,j = 0,1) molecules in 
collision with ^He atoms at ultra-low energies. They 
found that vibrational quenching of CO{v ~ 2, j = 0, 1) 
in collisions with ''He is dominated by the v = 2 ^ v = I 
transition. The cross sections for the w = 2 — s> w = 
transition were found to be about four orders of mag- 
nitude smaller than the single quantum transitions for 
both j = and j = 1 initial rotational levels. 

Ultracold vibrational relaxation of a number of other 
molecules in collisions with He atoms has been reported 
by a number of other investigators in recent years. 
Stoecklin, Voronin, and Rayez [38] reported the vibra- 
tional relaxation of F2 in collisions with '^He atoms. 
In this study, they constructed the PES of He— F2 us- 
ing ab initio points obtained by high-level molecular 
electronic-structure calculations, and reported the cross 
sections for elastic scattering and inelastic relaxation of 
F2{v = 0,1, j — 0) for collision energies in the range 
10~^ — 2000 cni^'^. A similar study has been reported 
by the same authors for the ■^He -I- HF(w = 0, 1, j = 0, 1) 
system (39| . The vibrational quenching cross sections 
were found to be very small compared to pure rota- 
tional quenching, in agreement with the results for the 
He— CO system. This is due to the weak dependence of 
the He— HF PES on the HF internuclear distance and the 
strong anisotropy of the interaction potential. 

Bodo and Gianturco [1^ presented a comparative 
study of vibrational relaxation of CO{v = 1,2, j — 0), 
HF(i; = l,2,j = 0) and LiH(w = l,2,j = 0) in colli- 
sions with "^He and ^He atoms in the Wigner regime. 
The quenching rate coefficients were found to depend 



strongly on the collision partner. They reported rate 
coefficients in the range 10~^^ — lO"'^^ cm'^ s""' for CO, 
10-1^-10-15 cm-'' s-i for HF, and lO-i-'-lO-" cm^ s'^ 
for LiH. The differences were attributed to the features of 
the intcrmolecular forces between the diatomic molecules 
and the He atoms. The interaction potential of the 
He— CO system is almost isotropic and is characterized 
by small vibrational couplings elements. The He— HF 
system is more anisotropic and the couplings between 
vibrational states are more significant. The interaction 
potential of the He— LiH system is very anisotropic and 
it exhibits strong vibrational couplings. 

There is considerable ongoing experimental interest 
in cooling and trapping NH [11 [il and OH [H, [13] 
molecules using the buffer gas cooling and Stark decel- 
eration methods described in the chapters by Doyle and 
by Meijer. Krems et al. [i^ and Cybulski et al. 
reported cross sections and rate coefficients for elastic 
scattering and Zeeman relaxation in '^He— NH collisions 
from ultralow energies to 10 cm~i. The calculations 
were performed using the rigid rotor approximation and 
an accurate He— NH PES. It was demonstrated that the 
elastic scattering of NH molecules with He atoms in weak 
magnetic fields is at least five orders of magnitude faster 
than the Zeeman relaxation, which suggests that the NH 
molecule is a good candidate for buffer-gas cooling. In 
a related study Gonzalez-Sanchez et al. [1^ examined 
rotational relaxation and spin-flipping in collisions of 
OH with He atoms at ultralow energies. They found 
that the rotational relaxation processes dominate the 
elastic process as the collision energy is decreased to zero. 

While theoretical prediction of rate coefficients for 
vibrational and rotational relaxation in a number of 
atom - diatom systems has been made, comparable 
experimental results are not available for majority of 
these systems. The first measurements of vibrational 
relaxation of molecules at temperatures below 1 K were 
reported by Weinstein et al. [4g. In their study, CaH 
molecules slowed down by elastic collisions with '^He 
buffer gas atoms were trapped in an inhomogeneous 
magnetic field. An upper bound of the rate coefficients 
for spin-flipping transitions in CaH as well as vibrational 
relaxation of CaH molecules in the v — 1 vibrational 
level in collisions with ■^He atoms were estimated at 
a t emp erature of about 500 mK. Balakrishnan et 
al. [43] presented a theoretical analysis of the vibra- 
tional relaxation of CaH in collisions with ^^He atoms 
based on quantum close-coupling calculations and an 
ab initio PES for the He— CaH system developed by 
Groenenboom and Balakrishnan [50] ■ In a related 
study, Krems et al. [5l| reported cross sections for 
spin-flipping transitions in CaH induced by collisions 
with '^He and obtained results in close agreement with 
the experimentally derived values of Weinstein et al [i^ . 
Krems et al. demonstrated that at low energies, spin- 
flipping transitions in the = rotational level of 
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molecules induced by structureless atoms occur through 
coupling to the rotationally excited iV > levels and 
that the corresponding rate coefficients are determined 
by the spin-rotation interaction with the transiently 
rotationally excited molecule. 

Table U provides a compilation of zero-temperature 
quenching rate coefficients for vibrational and rotational 
relaxation in a number of atom - diatom systems. 



system initial {v,j) A:r=o (cm^s ^) Ref. 
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TABLE L Zero-temperature inelastic rate coefficients for dif- 
ferent atom - molecule systems. 



2. Shape resonances in molecular collisions 

At energies above the onset of the s-wave regime, cross 
sections will be dominated by contributions from non- 
zero angular momentum partial waves. If the interac- 
tion potential includes an attractive part, the effective 
potentials for non-zero angular momentum partial waves 
may possess centrifugal barriers that introduce shape res- 
onances in the collision energy dependence of the cross 
section. This is illustrated in Fig|4] for the vibrational 
relaxation of C0(?; = l,j = 0) in collisions with "^He 
atoms. 

The sharp features in the energy dependence of the 
cross section for energies between 0.1 and 10.0 cm^^ 
arises from shape resonances supported by the van der 
Waals interaction potential between He and the CO 
molecule. As shown in FiglSJ when integrated over the ve- 
locity distribution of the colliding species the shape reso- 
nances lead to significant enhancement of the vibrational 




Energy (cm ') 

FIG. 4: Cross section for the quenching of the v = l,j = 
level of CO in collisions with *He as a function of the incident 
kinetic energy. Reproduced with permission from Balakrish- 
nan et al. [2311. 



relaxation rate coefficient for temperatures between O.I 
and 10.0 K. Similar results have been found for vibra- 
tional relaxation of CO Oz^, and CaH in 
collisions with ■^He atoms. The sharp features in the en- 
ergy dependence of the vibrational relaxation cross sec- 
tions for the Ar -I- D2 sytem shown in Fig[2] also arises 
from shape resonances supported by the Ar— D2 van der 
Waals potential. The effect is generally more pronounced 
for systems composed of heavier diatomic molecules and 
interaction potentials with deep van der Waals wells for 
which the density of states will be much higher leading 
to rich resonance structures in the cross sections. 



3. Feshbach resonances in molecular collisions 

Feshbach resonances occur in multichannel scattering 
in which an unbound (continuum) channel is coupled to 
a bound state of another channel. If the energy of the 
interacting system in the unbound channel lies close to 
that of the bound state and the coupling between the two 
channels is strong the cross section may change dramat- 
ically in the vicinity of the resonance. In the Feshbach 
resonance method for producing ultracold molecules, an 
external magnetic field is used to tune the energy of the 
bound pair to that of the the separated atoms. In atom 
- diatom systems, the bound state may correspond to a 
quasibound state of the atom - diatom van der Waals 
complex. Channel potentials corresponding to different 
initial vibrational and rotational levels of the diatom may 
induce Feshbach resonances. For the He— CO system 
Feshbach resonances were found to occur near channel 
thresholds corresponding to the j = 1 rotational level 



7 




T(K) 

FIG. 5: Rate coefficients for tfie quencliing of CO(« — l,j — 
0, 1) by coilisions with ^He as functions of tlie temperature. 
Reproduced witli permission from Balakrislman et al. [23| . 



in the v = and v = 1 vibrational levels. FiglS] shows 
the Feshbach resonance in the elastic scattering cross sec- 
tions in the v = 1, j ~ channel in the vicinity of the 

V = l,j = 1 level. The presence of the Feshbach reso- 
nance close to the opening of the j = 1 level has a dra- 
matic effect on the vibrational quenching cross sections 
from the v = l,j = 1 level of the CO molecule. Since 
the Feshbach resonance occurs so close to the threshold 
of the V = l,j ~ 1 channel, its effect on scattering in 
the ?; = l,j = 1 level is similar to that of the zero- 
energy resonance discussed previously for the Ar + D2 
system. This is illustrated in Fig. [S] (see also Table |l| 
where we compare the rate coefficients for vibrational 
relaxation from the v = l,j ^ and v — l,j — 1 lev- 
els of the CO molecule. The zero-temperature limiting 
value of the quenching rate coefficient of the v = 1, j = 1 
level is about two orders of magnitude larger than for the 

V = 1, j = level. Similar Feshbach resonances have also 
been shown to occur in the vibrational and rotational 
predissociation of He— H2 van der Waals complexes [20j . 
Forrey et al. [l^l has successfully used the effective range 
theory to predict the predissociation lifetimes of these 
resonances. 

The Feshbach resonances can be used as a very sensi- 
tive probe for the interaction potential and also to selec- 
tively break or make bonds in chemical reactions. The 
coupling between the bound and unbound states can be 
modified by applying an external electric or magnetic 
field and this provides an important mechanism for creat- 
ing or eliminating Feshbach resonances and thereby con- 
trolling the coUisional outcome. Krems have shown that 
weakly bound van der Waals complexes can be dissoci- 
ated by tuning a Feshbach resonance using an external 
magnetic field [5^. In this case the dissociation occurs 
through coupling between Zeeman levels of the bound 



and unbound channels and the magnitude of the cou- 
pling is varied by changing the external magnetic field. 

260 I ' ^ ' . 1 




200 I ' ' ' ' ^ ' 

3.8080 3.8090 3.8100 3.8110 

Energy (cm ) 

FIG. 6: Feshbach resonance in the elastic scattering cross 
section of CO{v = 1, j — 0) by *He atoms. The resonance 
occurs just below the opening of the v = 1, j — 1 level shown 
by the vertical line. The energy is relative to the v = 1, j = 
level of the CO molecule. Reproduced with permission from 
Balakrishnan et al. [23l |. 



B. Quasi- resonant transitions 

While the properties of cold and ultracold collisions 
are quite different from scattering at thermal energies 
and quantum effects dominate at low temperatures, a 
remarkable correlation between classical and quantum 
dynamics has been discovered in the relaxation of ro- 
vibrationally excited diatomic molecules. Experiments 
performed nearly two decades ago [H, [13] showed that 
collisions of rotationally excited diatomic molecules with 
atoms may result in very efficient internal energy trans- 
fer between specific rotational and vibrational degrees 
of freedom. The energy transfer becomes highly effi- 
cient when the collision time is longer than the rota- 
tional period of the molecule. This effect has since been 
termed "quasi-resonant rotation- vibration energy trans- 
fer". The experimental results revealed that the quasi- 
resonant (QR) transitions satisfy the propensity rule 
Aj = —AAv or A7' = —2 Aw where Av = Vf — vt and 
= if ~ ji [H: 13 ■ This inelastic channel dominates 
over all other ro- vibrational transitions. The QR transfer 
is generally insensitive to details of the interaction poten- 
tial. Rather, the QR process involves conservation of the 
action, I = riyV + njj, where Uy and rij are small inte- 
gers. Forrey et al. [2l| found that the QR transitions also 
occur in cold and ultracold collisions of rotationally ex- 
cited diatomic molecules with atoms and that the process 
is largely insensitive to the details of the interaction po- 
tential even in the ultracold regime. The Aj = —2Av QR 
transition in He + H2 collisions [55|] is illustrated in Fig. [7] 
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where the zero-temperature vibrational and rotational 
transition rate coefficients for different initial vibrational 
levels of the H2 molecule are plotted as functions of the 
initial rotational level. For initial rotational levels greater 
than 12 the QR transition becomes the dominant energy 
transfer mechanism compared to pure rotational quench- 
ing. The gap at j = 22 occurs because the Aj — —2Av 
transition is energetically not accessible for this initial 
state at zero temperature. Forrey et al. [2lj found that 
the QR process is even more dominant at low temper- 
atures than at thermal energies. Remarkably, classical 
trajectory calculations [2lj were successful in correctly 
predicting the correlation between Aj and Aw at ultra- 
cold temperatures even though the changes in v and j 
were fractional. Extensive studies of QR energy transfer 
in cold and ultracold temperatures have been reported by 
Forrey et al. [56l - l59| . Ruiz and Heller [g^] have recently 
published a review paper providing a detailed analysis of 
QR phenomenon u sing semi-classical techniques. McCaf- 
fery and coworkers |6ll - [63| have also reported a number of 
quasi-classical trajectory calculations of the QR process 
in thermal energy collisions and successfully interpreted 
a large body of experimental data based on the QR phe- 
nomenon and simple parametric models. 




i 

FIG. 7: Zero-temperature rate coefficients for *He + H2{v,j) 
coUisions as functions of the initial vibrational and rotational 
quantum numbers. Reproduced with permission from Forrey 
et al. [si]. 



C. Atom - molecular ion collisions 

Dynamics of ionic systems are different from collisions 
of neutral species. The short-range part of the interaction 
potential for ionic systems is usually more anisotropic 
and the long-range part has an attractive component 



which is determined by the polarizability of the atom and 
it vanishes as 1 /R* where R is the atom - molecule center- 
of-mass distance Due to the strong polarizability 

term, the interaction potential for ionic systems extends 
to longer range compared to neutral atom - molecule sys- 
tems. Therefore, it is important to understand the effect 
of both the short and long range part of the interaction 
potential on the scattering dynamics. For this purpose, 
several studies have focused on ultracold collisions be- 
tween molecular ions and neutral atoms as well as neutral 
molecules and atomic ions. 

Bodo et al. ^64ij investigated rotational quenching 
in NeJ + Ne and HeJ + He collisions at ultra-low 
energies. They found that the Wigner regime begins at 
a collision energy of 10~^ cm^^ for the He system and 
10~^ cm"-*^ for the Ne system. In general, the s-wave 
Wigner regime was found to occur at lower energies for 
ionic systems compared to neutral species. For example, 
in He -I- H2 [Hi and He -I- O2 [11] collisions the Wigner 
regime begins at collision energies of about 10^^ cm^^. 
The differences are attributed to the long range of 
the ion - neutral interaction potential which enhances 
contributions from higher partial waves. The differences 
between the He and Ne systems can be attributed to the 
mass difference and to the strength of the long range 
interaction potentials. For the heavier Ne system the 
long range part is more attractive, which increases the 
contribution of higher-order partial waves at ultra-low 
energies. The magnitude of the zero-energy rate coeffi- 
cients for rotational quenching in these molecular ions 
is on the order of 10^^ cm'^ s^^, which is considerably 
larger than for collisions involving neutral species. 

Stoecklin et al. |63| have performed a comparative 
study of collisions of NJ(d = l,j = 0) molecular ions 
with neutral '^He or ''He atoms, and the scattering of 
neutral N2(f = 1, j = 0) molecules in collisions with ■^He 
or "^He atoms. The vibrational quenching cross sections 
for these systems are presented in Fig. [H] While the 
behavior of the quenching cross sections of neutral N2 
molecules with ■^He and ^He atoms was found to be 
similar, they observed a striking difference between the 
quenching cross sections of NJ in collisions with '^He 
and *He atoms. In particular, the resonance positions 
in the quenching cross sections were significantly shifted 
and the number of resonances was different. The NJ 
-|- ■^He system has a shape resonance at a collision 
energy of lO^'^ cm^'^ which is absent for the Nj -t- ''He 
system. Furthermore, the zero-energy quenching rate 
coefficient is an order of magnitude larger for collisions 
of N^(i> — l,j — 0) with '^He than with *He, whereas 
it is comparable for collisions of N2(i; = l,j = 0) with 
both -^He and ''He. The differences may be due to a 
virtual state in the N^(i; = l,j = 0) -I- ^He collision 
system. More recently, Guillon et al. [gI] studied the 
effect of spin-rotation interaction on vibrational and 
rotational quenching for the He— N^ system. They 
found that the vibrational quenching is not modified by 
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FIG. 8: Vibrational quenching cross sections for Nj(ii = 
l,i = 0) + ^He and NJ(« = 1, j = 0) + "He (upper panel) 
and N2(t; = 1, j = 0) + ^He and N2(i; = 1, j = 0) + "He 
(lower panel) systems. Reproduced with permission from 
Stoecklin et al. 16511. 
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TABLE H: Zero-temperature inelastic rate coefficients for dif- 
ferent ion atom - molecule systems. 

isotopic analogues [13, [sl] . Most experimental studies of 
small molecular systems focussed on reactions at thermal 
or elevated collision energies though some recent mea- 
surements have been extended to temperatures as low as 
10 K for some astrophysically relevant systems [69]. Due 
to the possibility of achieving coherent chemistry there 
has been substantial interest in understanding the behav- 
ior of chemical reactions at cold and ultracold tempera- 
tures. At these temperatures perturbations introduced 
by external electric and magnetic fields are significant 
compared to the collision energies involved and external 
fields can be employed to control and manipulate the 
reaction outcome. Over the last seven years a number 
of studies of ultracold atom - diatom chemical reactions 
have been reported both for reactions with and without 
an energy barrier. Chemical reactions at low temper- 
atures often behave quite differently from reactions at 
elevated temperatures. In particular, the weak van der 
Waals interaction potential which does not play any sig- 
nificant role at high temperatures may have a dramatic 
effect on the outcome of reactions at low temperatures. 



the spin-rotation coupling, while rotational transitions 
are sensitive to the fine-structure interactions. 

Table |TT] provides a compilation of zero-temperature 
quenching rate coefficients for vibrational and rotational 
relaxation in several ion atom - molecule systems. 



III. CHEMICAL REACTIONS AT ULTRACOLD 
TEMPERATURES 

The last three decades have seen impressive progress 
in the experimental and theoretical descriptions of chem- 
ical reactions between atoms and small molecular sys- 
tems. While fully state-resolved experiments have been 
performed for a large number of collision systems, accu- 
rate quantum dynamics calculations have been restricted 
to systems involving light atoms such as H -I- H2, F -t- 
H2, CI -I- H2, C -I- II2, N -I- II2, O -f H2 and some of their 
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A. Tunneling dominated reactions 

In the course of the last 10 — 15 years there has 
been much interest in understanding the role of res- 
onances in chemical reactions that involve an energy 
barrier for which the reactivity is primarily driven by 
tunneling at low temperatures. Recent studies on F 
+ H2/HD/D2 P IIOH, CI + HD 

HCl/DCl m, Li + RF / LiF + H ^ 

H2 [8l,[83, and F + HCl/DCl 'S^l reactions have demon- 
strated that decay of quasibound states of the van der 
Waals interaction potential in the entrance and exit chan- 
nels may give rise to narrow Feshbach resonances in the 
cross sections. 

The reactions of F with H2 and HD have been the sub- 
ject of numerous quantum scattering calculations over 
the last two decades. The two reactions have emerged 
as benchmark systems for experimental and theoreti- 
cal studies of resonances in chemical reactions. A de- 
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tailed analysis of the low energy resonance features in 
the F + H2 reaction was presented by Castillo et al. in 
1996 [70]. They demonstrated that several resonances 
that appear in the energy dependence of the cumulative 
reaction probability for the F + H2 reaction arise due to 
the van der Waals interaction potential in the product 
HF + H channel. In particular, the origin of the reso- 
nances has been attributed to the van der Waals potential 
associated with the HF(ti' = 3, / = — 3) channels. A 
large number of experimental and theoretical papers has 
appeared which examined various aspects of these reso- 
nance features [TTI. [70l - l72l . [75l - l80| . Among these studies, 
the work of Takayanagi and Kurosaki [7l[ deserves spe- 
cial attention. They showed that for F -I- H2, F + HD 
and F -I- D2 reactions reactive scattering resonances oc- 
cur due to the decay of quasibound states of the van der 
Waals potential in the entrance channel of the reaction. 
These Feshbach resonances are associated with the decay 
of quasibound states of adiabatic potentials correspond- 
ing to F- • • H2(v = 0, J = 0, 1), F- • • HD(« - 0,j - 0-2), 
and F---D2(u — 0,j — — 2) complexes, obtained 
by diagonalizing the J = Hamiltonian in a basis set 
of the asymptotic ro-vibrational states of the reactant 
molecules. 




Kinetic energy (eV) 



FIG. 9: J = cross sections for F + H2(w = 0,j = 0) 
IIF(u') + H reaction for v' = — 2 as functions of the incident 
kinetic energy. Reproduced with permission from Balakrish- 
nan and Dalgarno [lH . 



Reactions at zero temperature 



Balakrishnan and Dalgarno [11[ showed that quasi- 
bound states of the F -t- II2 van der Waals complex have 
a dramatic effect on the reactivity in the Wigner thresh- 
old regime. They found that the F -|- H2(v = 0,j = 0) 
reaction has a rate coefficient of 1.25 x 10~^^ cm^ s~^ in 
the zero-temperature limit. Their calculation was based 
on the widely used PES of Stark and Werner [s^ for 
the F -I- H2 system. The relatively large value of the 
zero-temperature rate coefficient is due to the presence 
of a small narrow energy barrier for the reaction so that 
tunneling of the H atom is efficient. The vibrational 
level resolved cross sections for the F -f H2(w = 0, = 0) 
reaction are shown in Fig. IH] for the total angular 
momentum quantum number J = 0. The v' = 2 level is 
the dominant channel at low energies in agreement with 
the behavior at higher energies. 

Subsequent calculations showed that at low energies 
the F -I- HD reaction is dominated by the formation 
of the HF product with an HF/DF branching ratio of 
about 5.5 [74|. The formation of the DF product is 
suppressed because tunneling of the heavier D atom 
is less efficient. In earlier quantum calculations, Baer 
and coworkers [sO, [sH reported HF/DF branching 
ratios ranging from 1.5 at 450 K to 6.0 at 100 K. Their 
low temperature value is in close agreement with the 
zero-temperature limit obtained by Balakrishnan and 
Dalgarno [t^] . Fig. [10] compares the J — cross sections 
for the F + YL2{v = 0,j = 0) and F + HD(t; = 0,.? = 0) 
reactions over a kinetic energy range of 10"^ to 1.0 eV. 



In the Wigner threshold regime, the reactivity of the 
F -|- H2 system is an order of magnitude greater than 
that of the F + HD reaction. A rigorous analysis 
of the scattering resonances in the F -I- HD reaction 
was recently presented by De Fazio et al. [s^, who 
provided a detailed characterization of the resonances 
supported by the entrance and exit channels of the van 
der Waals potential and discussed the effect of higher 
total angular momenta on the position and lifetime of 
the resonances. Stereodynamical aspects of the F -|- 
H2 collision and the effect of polarization of the H2 
molecule on the outcome of the reaction at low energies 
were recently explored by Aldegunde et al. (tJ- They 
argued that a reactant polarization scheme can be ex- 
ploited to control state-to-state dynamics of the reaction. 

To explore the role of tunneling in chemical reactions 
at cold and ultracold temperatures Bodo, Gianturco, and 
Dalgarno [t^ investigated the dynamics of the F -f D2 
system at low and ultra-low energies. They found that 
compared to the F -I- H2 reaction, the reactivity of F 
-|- D2 is significantly suppressed in the Wigner regime 
with an HF/DF ratio of about 100. This is illustrated in 
Fig. [TT] where the J = cumulative reaction probability 
and cross sections for the F -I- H2 and F -I- D2 reactions 
are plotted as functions of the incident kinetic energy 
from the Wigner limit to 0.01 eV. 

The dramatic suppression of the F -I- D2 reaction in the 
Wigner regime cannot be explained based on tunneling 
alone. A closer examination of the reaction probabilities 
for the F -I- H2 and F -|- D2 reaction revealed that an un- 
usual enhancement of the reactivity occurs for the F -|- 
H2 reaction at about 3 x 10~^ eV [7a| (see the upper panel 
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FIG. 10: Comparison of J = cross sections for F + HD(« — 
0,i = 0) -5> HF + D (solid curve), F + HD(u = 0, j = 0) 
^> DF + H (sliort-daslied curve) and F + U2{v = 0,j = 0) 
HF + H (long-dashed curve) reactions as functions of the 
incident kinetic energy. Reproduced with permission from 
Balakrishnan and Dalgarno 



of Fig. [TT|) . This feature was attributed to the presence 
of a virtual state. The enhancement of the hmiting value 
of the rate coefficient for the F + H2 reaction was also as- 
cribed to the virtual state. Evidence for the virtual state 
is the presence of a Ramsauer-Townsend minimum in the 
elastic cross section at an energy of about 3 x 10~^ eV 
and a negative value of the real part of the scattering 
length for the F -f- H2 reaction j75| . 

To explore how isotope substitution modifies reac- 
tivity at low temperatures Bodo et al. [75] artificially 
varied the mass of the hydrogen atoms in the calculation 
for the F -f- H2 reaction from 0.5 to 1.5 amu. As 
illustrated in Fig. [TH for an H atom mass of 1.12 amu 
the virtual state induces a zero energy resonance at 
which the real part of the scattering length diverged 
to infinity. For the same value of the H atom mass, 
the zero-temperature rate coefficient of the reaction 
attains a value of 1.0 x 10~^ cm^ s^^ which is about 
three orders of magnitude larger than that of the F + 
H2 reaction in the Wigner limit. The variation of the 
scattering length of the F -I- H2 system as a function 
of the mass of the pseudo-hydrogen atom is similar 
to the variation of scattering length as a function of 
the magnetic field in the vicinity of a Feshbach resonance. 

Another example of a tunneling dominated reaction is 
the Li -I- HF LiF + H reaction. At cold and ultracold 
temperatures the reaction occurs by tunneling of the rel- 
atively heavy fluorine atom. The LiH -|- F channel is 
energetically not accessible at low energies and tunnel- 
ing of the H atom is not involved in this reaction at low 
temperatures. Due to strong electric dipole forces ex- 
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FIG. 11: Comparison of the J = cumulative reaction prob- 
abilities (upper panel) and cross sections (lower sub-panel) of 
F -I- H2 and F -I- D2 reactions as functions of the incident 
kinetic energy. Reproduced with permission from Bodo et 
al. M. 



erted by the HF molecule, the van der Waals interaction 
potential of the LiHF system is deeper compared to the 
F -I- H2 system. The Li- • ■ HF van der Waals potential 
weU is about 0.24 eV (1936.0 cm~i) and the H- • • LiF 
potential well is about 0.07 eV (565.0 cm^^). Quantum 
scattering calculations for Li + HF and LiF + H reac- 
tions by Week and Balakrishnan jMi^] with the PES of 
Aguado et al. [92| have shown that for incident energies 
below 10"'^ eV the reaction cross section exhibits a large 
number of resonances. The energy dependence of the 
J = cross section for Li -I- HF(w — 0,j — 0) collisions is 
shown in Fig. [131 Detailed bound state calculations of the 
LiHF van der Waals complexes revealed that for the Li 
+ HF(i; = 0, j = 0) reaction the resonances correspond 
to the decay of Li- • • HF(i; = 0, j = 1 — 4) van der Waals 
complexes. Calculations with vibrationally excited HF 
molecules showed that the reaction becomes about 600 
times more efficient in the Wigner regime when HF is 
excited to the v = 1 vibrational level. As seen in Fig. [T31 
a unique feature of the Li -|- HF(v — Q, j ~ 0) reac- 
tion is the presence of a strong peak at 5.0 x 10^'' eV 
at which the reaction cross section is about six orders of 
magnitude larger than the background cross section. The 
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FIG. 12: Real part of the scattering length (upper panel), 
zero-temperature limiting values of the rate coefficient (mid- 
dle panel), and positions of quasibound states (lower panel) 
of the F- • • H2 complex as functions of the mass of a pseudo 
hy drog en atom. Reproduced with permission from Bodo et 
al. IM. 



results for both Li + HF |8J and LiF + H [85| reactions 
with thermal and non-thermal vibrational excitation sug- 
gest that heavy-atom tunnehng may play an important 
role in chemical reactions at cold and ultracold temper- 
atures. An experimental study [l^] of an organic ring 
expansion reaction at 8 K has shown that the reaction 
occurs almost exclusively by carbon tunneling. The tun- 
neling contribution was found to be orders of magnitude 
greater than over the barrier contribution. In a recent 
work, Tscherbul and Krems [9^ have explored the Li -|- 
HF and LiF -I- H reactions in the presence of an external 
electric field. They have shown that, for temperatures 
below 1 K, the reaction probability can be significantly 
influenced by electric fields. 



2. 



Feshbach resonances in reactive scattering 



Van der Waals complexes formed during collisions can 
either undergo vibrational predissociation or vibrational 
prereaction leading to sharp features in the energy de- 
pendence of the cross sections. The term "prereaction" 
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FIG. 13: Cross sections for LiF formation (solid curve) and 
non-reactive scattering (dashed curve) in Li -|- HF{v,j = 0) 
collisions as functions of the incident kinetic energy: results 
for V — (lower panel); results for u = 1 (upper panel). Re- 
produced with permission from Week and Balakrishnan [83 |. 



refers to the process in which a rotationally or vibra- 
tionally excited van der Waals complex decays through 
chemical reaction rather than by rotational or vibrational 
predissociation. For reactions with energy barriers the 
chemical reaction pathway may involve tunneling. The 
reactions of CI with H2 and HD are dominated by tun- 
neling at low temperatures (sil - lssj . Compared to the F 
-I- H2 and F -I- HD reactions, the energy barrier for the 
CI + H2 reaction is much larger and the reactivity at 
low energies is significantly suppressed. Nearly a decade 
ago, Skouteris et al. [8l| showed that the van der Waals 
interaction potential between CI and HD determines the 
reaction outcome, despite the fact that the depth of the 
van der Waals interaction potential for the CI- • • HD sys- 
tem is less than one-tenth of the height of the reaction 
barrier. Quantum scattering calculations of the CI -I- 
HD reaction on PESs without the van der Waals interac- 
tion potential predict nearly equal probabilities for HCl 
and DCl products 
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However, if the potential surface 
includes the van der Waals interaction a strong prefer- 
ence for the DCl product occurs at thermal energies, in 
agreement with experimental results. The effect of these 
weakly bound states on the reactivity at cold and ultra- 
cold temperatures has recently been explored by Balakr- 
ishnan [82i] ■ 

The cross sections for HCl and DCl formation and non- 
reactive rovibrational transitions in CI + HD(u = l,j = 
0) collisions for total angular momentum quantum num- 
ber J ~ are shown in Fig. [14] as functions of the to- 
tal energy [82i|. The sharp features in the cross sections 
correspond to Feshbach resonances arising from the de- 
cay of quasibound van der Waals complexes formed in 
the entrance channel of the reaction. The quasibound 
states can be identified by examining bound states of 
adiabatic potentials correlating with the w = l,j = 
and V = l,j = 1 levels of the HD molecule. They are 
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FIG. 14: Reactive and non-reactive scattering cross sections 
for CI + HD(« = l,j = 0) reaction as functions of the incident 
kinetic energy. Reproduced with permission from Balakrish- 
nan 18211. 



displayed in Fig. [15] as functions of the atom - molecule 
separation. The adiabatic potential curves are computed 
by diagonalizing the diabatic potential energy matrix ob- 
tained in a basis set of rovibrational levels of the HD 
molecule at each value of the atom - molecule separa- 
tion. The Feshbach resonances labeled B, C, D, and E in 
Fig. [14] result from the decay of the corresponding qua- 
sibound states shown in Fig. 1151 




The cross sections in Fig. [14] do not show a peak 
corresponding to the metastable state A. This is because 
the state A is too deeply bound and it is not accessible 
through scattering in the v = 1, j = channel. Fig. [T3] 
also shows that the quasibound states preferentially 
undergo prereaction than predissociation. This is 
due to the enhanced coupling of the resonance states 
with the reactive channels compared to those with the 
non-reactive channels. The wavefunctions of quasibound 
states B and E are shown in Fig. [JH] as functions of the 
atom - molecule separation. Although the wavefunction 
of the weakly bound state E extends far beyond the 
transition state region of the reaction, it preferentially 
undergoes prereaction compared to predissociation. 
Thus, regions of the interaction potential far away from 
the transition state region may have a significant effect 
on reactivity especially when part of the wavefunction in 
that region is sampled by a resonance state which is cou- 
pled to the reactive channel. Fig. [T3] also demonstrates 
that since the reaction is dominated by tunneling at low 
temperatures the formation of the DCl -I- H product 
which involves the tunneling of the D atom is severely 
suppressed in the threshold reg ime. This is more clearly 



illustrated in Fig. 4 of Ref. (82| where the reaction cross 



section is plotted as a function of the incident kinetic 
energy. 
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FIG. 16: Adiabatic potential and the wavefunctions of the 
quasibound levels B and E shown in Fig llSl as functions of the 
atom - molecule separation. Amplitudes of the wavefunctions 
have been reduced by a factor of 10 for the convenience of plot- 
ting. Reproduced with permission from Balakrishnan [82l ]. 



FIG. 15: Adiabatic potential energy curves of the CI + HD 
system correlating with the IID(ii — l,j = 0) and HD(u = 
1, J = 1) levels as functions of the atom - molecule separation. 
Quasibound levels responsible for the resonances in Fig |14l are 
labeled by B, C, D, and E. Reproduced with permission from 
Balakrishnan [821. 



Table IIIII provides a compilation of zero-temperature 
quenching rate coefficients for a number of atom - diatom 
chemical reactions which are dominated by tunneling at 
low energies. 
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TABLE IH: Zero-temperature quenching rate coefficients for 
tunneling dominated reactions. The final arrangement has 
also been specified. 



B. Barrierless reactions 

1. Collision systems of three alkali metal atoms 

As discussed by Hutson and Soldan in two recent re- 
views @, [3], progress on the production of ultracold 
molecules and the creation of molecular Bose-Einstein 
condensates of alkah-metal systems have motivated the- 
oretical studies on ultracold atom - dimer alkali-metal 
collisions. Here, we give an overview of recent cal- 
culations for spin-polarized triatomic alkali-metal sys- 
tems, Li + Lis ^3^, Na + Na,2 and K + 

K2 [lOOj |. These results have all been obtained using 
the reacti ve sc attering code written by Launay and Le 
Dourneuf |lOl| . based on a ti me- i ndep endent quantum 
formalism. Refs. [il [H [93, [lo3| describe the details 
of the PESs and dynamics calculations. The quantum 
dynamics studies of Refs. [98, 99] have been per formed 
using the Naa PES reported by Higgins et al. |l02| , while 
the results of Ref. [93| have been obt aine d using the Lis 
PES calculated by Colavecchia et al. [Toaj . Another PES 
for Lis has been constructed by Brue et al. jl04| . 

Quantum scattering calculations show that ultracold 



reactions of alkali metal atoms with alkali metal dimers 
are much more efficient than tunneling driven processes. 
This can be explained based on the following consid- 
erations. First, the indistinguishability of the atoms 
may play a role. For a homonuclear triatomic system, 
the three possible arrangement channels are the same. 
Also the presence of identical two-body potentials in all 
three arrangement channels may enhance the depth of 
the three-body interaction potential. Consider the two- 
body term (additive term) of a homonuclear triatomic 
system at equilateral geometries. Since the distances be- 
tween the three identical atoms are the same, the same 
diatomic potential term is added three times to yield the 
two-body terms of the triatomic system. If the diatomic 
potential is deep or repulsive, the two-body term will be 
three times deeper or repulsive at equilateral geometries. 
In contrast, for triatomic systems with distinguishable 
atoms, the diatomic pairs are not the same. For example, 
for the Li -I- HF system, the three diatomic fragments are 
distinctly different. They have completely different elec- 
tronic structures and properties such as minima, equilib- 
rium distances, turning points, and the nature and range 
of the interaction. Thus, at an equilateral configuration 
for which the three diatomic distances are the same but 
the diatomic potential energies are different (one may 
be attractive while the other two may be repulsive), the 
overall two-body term can be weaker than the individual 
two-body interaction potential. This could lead to a tri- 
atomic interaction region that is less strong than that of 
a homonuclear triatomic system. 

Second, the topology of the PES plays a significant 
role. For all alkali-metal trimer systems, the minimum 
energy configurations arise at equil atera l and collinear 
geometries as shown by Soldan et al. (l05| . Furthermore, 
the surface is barrierless so that all atom - dimer alkali- 
metal collisional approaches are energetically possible. In 
contrast, most of the non- alkali systems discussed above 
are characterized by a collinear (or bend) atom - diatom 
approach, and dominated by a repulsive barrier in the 
triatomic transition-state region. The particular topol- 
ogy of the PES arises from the three-body interaction 
potential. Also, the couplings between different elec- 
tronic surfaces of the triatomic system give rise to coni- 
cal intersections, which create a repulsive barrier in the 
transition-state region. If the energy barrier is high and 
its width large, tunneling will be very inefficient leading 
to very small values for the reaction rate coefficients in 
the Wigner regime. Although reactivity can be enhanced 
when resonances are present, the background scattering 
in reactive cross sections is generally quite small. In 
contrast, for almost all alkali-metal trimer systems, the 
conical intersections arise at small interatomic distances 
where the PES is sufficiently repulsive that it plays no 
significant role at ultra-low energy and the reactivity is 
not influenced by such repulsive barriers. 

Third, the density of states of the system plays an 
important role. If the density of states is small as for 
light systems, the typical energy spacing is rather large. 
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For heavy systems such as alkali-metal systems, the 
density of states is very large and the narrow energy 
spacing can lead to very strong couplings, especially in 
the vicinity of avoided crossings. The strong couplings 
will lead to efficient energy transfer between different 
quantum states. This may also explain the relatively 
weak dependence of the vibrational quenching rate coef- 
ficients on the initial vibrational state of the molecule in 
atom - dimer alkali-metal collisions. 

Fig. [T7] provides a comparison between elastic and 
que nchi ng rate coefficients for "^^K -I- ^^K2{v = l,j = 
0) (lOOj l collisions at energies ranging from 10~^ K to 
10~^ K. The results include contributions from the total 
angular momentum quantum numbers J = — 5. For 
low vibrational states of the K2 dimer, quenching pro- 
cesses are more efficient than elastic scattering at ultra- 
cold energies. For example, at a temperature of 10~^ K, 
the quenching rate coefficient is about 10"""^° cm^ s~^ 
compared to 10^^"^ cm'^ s^^ for elastic scattering. The 
quenching processes lead to coUisional relaxation of the 
molecules to lower ro-vibrational states resulting in trap 
loss. 

Similar results have been obtained for Na -|- 
Na2 iaa, i99i and Li + U2 [H, IH, ^ collisions. The 
zero-energy quenching rate coefficients for these systems 
are found to be on the order of 10~^^ — 10~^° cm^ s~^. 
The quenching processes are more efficient than the 
elastic collisions at ultracold temperatures. Thus, 
quenching of vibrationally excited alkali-metal dimers in 
collisions with alkali-metal atoms occurs with significant 
rates at ultracold temperatures. The typical magnitude 
of these rate coefficients from the scattering calculations 
is in reasonable agreement with experi menta l results. 
In a recent experiment, Staanum et al. jl06^ reported 
a value of 9.8 x 10~^^ cm'^ s~^ for the relaxation of 
low-lying vibrational levels of Cs2 in collisions with Cs 
atoms at a temperature of 60 x 10~^ K. In a separate 
experiment, Zahzam et al. }10'<^ determined a quenching 
rate coefficient of 2.6 x 10~^^ cm'^ for the same sys- 
tem at a temperature of 40 x 10~^ K. Wynar et al. [l08j | 
estimated inelastic rate coefficients of 8 x 10~^ ^ cm'' 
for Rb -I- Rb2 collisions. Mukaiyama et al. jlQ9l | reported 
an inelastic rate coefficient of 5.5 x 10"-'^^ cm^ s~^ for 
collisions of Na atoms with Na2 molecules created b y the 
Feshbach resonance method while Syassen et al. [ll0j | 
obtained a value of 2 x 10~^° cm'' for collisions of 
Rb atoms with Rb2 molecules. 

At large atom - diatom separations, R, the interaction 
potential can be approximated by an effective potential, 
composed of a repulsive centrifugal term and the long 
range interaction potential. For alkali-metal trimer sys- 
tems, the atom - diatom long range potential is a van 
der Waals interaction potential and behaves as - Gq/R^ . 
The classical capture model (also known as the Langevin 
model ) ha s been shown to work quite well for these sys- 
tems |10G| | at certain energy regimes. The Langevin 
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FIG. 17: Elastic (upper panel) and quenching (lower panel) 
rate coefficients versus the collision energy for '^^K -|- ^^K2(v = 
l,j — 0) scattering. The rate coefficient from a capture model 
is also shown for the quenchin g pro cesses. Adapted with per- 
mission from Quemener et al. [lOOl ]. 



model is based on the assumption that if the total energy 
of the system can overcome the effective barrier, then the 
elastic probability is zero and the quenching probabil- 
ity is one. Otherwise, the elastic probability is one and 
the quenching probability is zero. Classically, the barrier 
prevents the atom and the molecule from accessing the 
region of strong coupling with the other channels. 

Results obtained using the Langevin model for the K 
-I- K2 system are shown in Fig. [17] along with the quan- 
tum results. As Fig.[T7]illustrates, three different regions 
can be distinguished for the K + K2 system. The first 
region, for collision energies below 10~^ K, corresponds 
to the Wigner regime where the threshold laws apply. 
In this regime, the quenching rate coefficient tends to a 
constant while the elastic component approaches zero as 
the square root of the collision energy. A quantum de- 
scription is required and classical models cannot describe 
the dynamics in this regime. The second region, above a 
collision energy of 10 ~^ K, corresponds to the Langevin 
regime. Here, the difference between the full quantum 
calculation and the classical model is within 10 %. Simi- 
lar results have been found for ^Li + ^Li2(f = 1,2, j — 0) 
and ^Li + ^Li2(u = l,2,j = 1) collisions by Cvitas et 
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al. [94| , who showed that in the Langevin regime the rate 
coefficients become independent of the vibrational level 
of the molecule. The third region is intermediate be- 
tween the Wigner and Langevin regimes, which typically 
corresponds to a quantum calculation with two or three 
partial waves. As shown in the lower panel of Fig. [17] 
the height of each J-resolved effective barrier depicted 
by vertical lines for J = 1 — 5 corresponds approximately 
to the maximum of each J-resolved quenching rate co- 
efficient. This indicates that the quenching process be- 
comes significant when the barrier height is energetically 
overcome. Whenever the elastic and quenching rate co- 
efficients have comparable values, the Langevin regime is 
reached, as shown in Fig. 1171 Since the quenching prob- 
ability is almost equal to one and the elastic probability 
is almost equal to zero, elastic and quenching transition 
matrix elements are both equal to one and the cross sec- 
tions and rate coefficients for the two processes become 
similar. 

The above analysis shows that at energies above the 
onset of the Wigner regime the Langevin model can cor- 
rectly describe the dynamics for barrierless atom - dimer 
alkali-metal collisions because of the strong inelastic cou- 
plings. However, the classical model is not suitable for 
tunneling dominated reactions. Accurate quantum dy- 
namics calculations are computationally demanding for 
heavier systems such as Cs -I- Cs2, Rb + Rb2, Rb + 
RbCs or Cs -I- RbCs and the classical model may be 
used to qualitatively describe the dynamics of these sys- 
tems at energies above the s-wave regime. The tempera- 
ture dependence of the rate coefficients predicted by the 
Langevin model is given by the simple formula in atomic 
units: 




where fi is the reduced mass of the atom - diatom colli- 
sional system and Ce the dominant atom - diatom long 
range coefficient. In Fig.[TH]we show the rate coefficients 
predicted by the Langevin model as functions of the tem- 
perature for different atom - diatom and diatom - di- 
atom collisions involving Rb or Cs atoms. For Rb— RbCs, 
Cs— RbCs, and RbCs— RbCs collisions we used the corre- 
sponding Cg coefficients calculated by Hudson et al. |lll| . 
In the absence of similar data for Cs -I- Cs2 and Rb -|- Rb2 
collisions, we approximated the atom - dimer Cq coeffi- 
cient by multiplying the corresponding atom - atom value 
by a factor of two. For Cs -I- Cs2 collisions we used the Cg 
coe fficie nt for Cs— Cs interaction calculated by Amiot et 
al. (ll2j | and for Rb -|- Rb2 collisions we adopted the Cg 
coefficient for Rb— Rb interaction reported by Derevianko 
et al. [ll3l |. As shown in Fig.[T8]the Langevin model pre- 
dicts rate coefficients within an order of magnitude of 
the experimental values for the different systems. The 
predicted results for the Cs -I- Cs2, Rb -I- RbCs, and Cs 
-f- RbCs collisions ag ree with t he c orresponding experi- 
mental data of Refs. [l06l llOTl . Illl| within the reported 
error bars. For Rb -I- Rb2 collisions the Langevin model 



predicts results in very close agreement with the experi- 
mental result of Wynar et al. jl08j |. Thus, the Langevin 
model appears to be valid for describing collisional prop- 
erties of these systems at /iK temperatures. The model 
also confirms that cold and ultracold collisions of alkali 
metal atoms and dimers are essentially characterized by 
the leading term in the long-range part of the interaction 
potential. 
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FIG. 18: Rate coefficients as functions of the temperature 
predicted by the Langevin model compared with experimen- 
tal data for different collisional processes involving Rb or Cs 
atoms. The curves correspond to the Langevin results and 
the symbols denote the experimental results. 



2. Role of PES in determining ultracold reactions 

Potential energy surfaces are the key ingredients that 
enter in quantum dynamics calculations. While the two- 
body terms are generally well known and accurate, the 
three-body terms are more difficult to compute with high 
precision since they are non-additive and involve correla- 
tions between the three atoms. As a consequence, quan- 
tum dynamics calculations may suffer from the quality 
and degree of accuracy of the three-body terms. The 
sensitivity of the ultracold collision cross sections to the 
details of the PES has been investigated for atom - dimer 
alkali-metal systems by Quemener et al. [o^ for Na -I- 
Na2(w = 1 - 3,j = 0) and by Cvitas et al. ^ for Li 
+ Li2{v = ~ 3,j = 0). In these studies, a linear scal- 
ing factor A has been included to tune the three-body 
interaction term in the PES and cross sections were cal- 
culated as a function of A. Other studies have com- 
pared the dynamics with and without the three-body 
terms for Na + Na2{v = l,j = 0) [98] and for Li + 
Li2(i; = - 10, j = 0) l96( collisions. 

The dependence of the total quenching cross sections 
on the three body term for Na -I- Na2(w = I — 3, j — 0) 
scattering [ooj at a collision energy of 10~^ K is shown in 
Fig. [m The contribution of each final vibrational level 
is also plotted. For the v = 1 vibrational level, the cross 
sections are very sensitive to the details of the three-body 
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FIG. 19: Dependence of the quenching cross sections on the 
three-body term at a collision energy of 10~^ K for Na + 
Na2(w = 1 — 3,j = 0). Reproduced with permission from 
Quemener et al. |9!|]. 



potential. A change of 1 % in A leads to a significant 
change of 75 % in the cross sections. At the minimum 
of the Na3 potential, a change of 1 % corresponds ap- 
proximately to 10 K. At present, ab initio calculations of 
PESs cannot be done with such accuracy. Thus it ap- 
pears that it is difficult to get an accuracy better than 
two orders of magnitude in the cross sections for Na -I- 
Na2(w — l,j = 0) collisions. However, for v — 2 and 
3, the total cross sections show a weaker dependence on 
the three-body term. The state-resolved cross sections 
do not show strong dependence on the three-body term, 
compared to collisions of molecules in u = 1. 

Fig.[20lshows the dependence of the state-to-state cross 
sections on the three body term for the reaction Na -|- 
Na2(w = 3,j = 0) Na Na2(w/ = 2,l,0,j/) at a col- 
lision energy of 10~^ K. The oscillations in the cross sec- 
tions, when A is modified, are due to Feshbach resonances 
which arise when a triatomic quasibound state (or a vir- 
tual state) crosses the energy threshold as the strength of 
the interaction potential is decreased (increased). When 
such a Feshbach resonance occurs, Cvitas et al. [93] have 
argued, strong resonant peaks appear in the cross sec- 
tions if inelastic couplings are weak. In contrast, weak 
oscillations of one order of magnitude at best appear 
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FIG. 20: Variation of state-to-state cross sections at a colli- 
sion energy of 10~^ K for Na -I- Na2(w = 3, j = 0) ^ Na -I- 
Na2(u/ = 2, 1,0,j7) as a function of the parameter A. See 
text for details. Reproduced with permission from Quemener 
et al. fii. 



when inelastic couplings are strong, as in alkali-metal 
trimer systems. A ge nera lization of this effect has been 
discussed by Hutson [ll4j . Larger modifications of the 
three-body term affect the cross sections significantly. 
This can also be seen in Fig.[2T]for Li -I- Li2 collisions [9^] 
in which the three-body term is excluded from the dy- 
namics calculation. The state-to-state cross sections ex- 
hibit a stronger dependence on the three body term. 
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3. Relaxation of vibrationally excited alkali metal dimers 

Ultracold diatomic molecules produced in pliotoas- 
sociation or Feshbach resonance methods are usually 
created in excited vibrational states. Theoretical studies 
involving highly vibrationally excited molecules are 
challenging due to the large number of energetically 
open reaction channels present in the quantum calcu- 
lation. This puts severe restriction on the calculations 
for vibrationally excited molecules at low temperatures. 
Ultracold quantum dynamics calculations for collisions 
of highly vibrationally excited molecules have been 
reported in 2007 by Quemener et al. [o^ for the Li + 
Li2 system. A three-atom problem is generally described 
by two different kinds of asymptotic channels. The 
first kind are the single-continuum states (SCSs). They 
correspond to configurations where the triatomic system 
dissociates asymptotically into an atom and a diatomic 
molecule. The second kind are the double-continuum 
states (DCSs). They correspond to configurations 
where the triatomic system dissociates asymptotically 
into three separated atoms. Since highly vibrationally 
excited molecular states lie close to and below the 
triatomic dissociation limit, they are also coupled with 
the DCSs which lie above the dissociation limit. As 
a consequence, DCSs have to be included in quantum 
simulations of atom - diatom systems involving highly 
vibrationally excited diatomic molecules [96|. This 
dramatically increases the size and complexity of the 
quantum dynamics problem. 
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FIG. 21: Dependence of the elastic and quenching rate coef- 
ficients on the vibrational excitation of the molecule for 'Li 
-I- 'Li2(t; = — 10, j = 0) scattering at a collision energy 
of 10~® K. The bold line corresponds to the full calculation 
and the thin line corresponds to the calculation without the 
three-body term of the PES. Adapted with permission from 
Quemener et al. JoE]. 



Figure [H] shows the dependence of the rate coeffi- 
cients for ^Li -I- ^Li2(f = — 10, j = 0) collisions on 



the vibrational quantum number at a collision energy 
of 10~^ K |96|. Quenching processes are more efficient 
than the elastic scattering for both high and low vibra- 
tional levels. Similar results have been found for ^Li + 
®Li2(u = — 9, j = 1) ^] system composed of fcrmionic 
atoms. Quenching rate coefficients show a slight de- 
crease when the molecule is in its highest vibrational 
state. This is because the overlap of the wavefunctions 
of highly excited diatomic molecules with low-lying 
vibrational levels is very small leading to small values 
for the interaction potential coupli ng m atrix elements 



between the initial and final states [115|. These results 
do not directly apply to ultracold molecules composed 
of fermionic atoms created near a Feshbach resonance. 
In these experiments, the quenching processes are 
suppressed because the atom - atom scattering length 
is tuned to large and positive values leading to efficient 
Pa uli b locking mechanism, as explained by Petrov et 
al. [ll6l |. In the theoretical study of Li -I- Li2 collisions, 
the atom - atom scattering length is small and negative 
and no suppression of the quenching processes is found 
for the molecule in the last vibrational state. Thus, the 
sign and magnitude of the Li— Li scattering length play a 
crucial role in the mechanism that suppresses quenching 
collisions. 

The quenching rates displayed in Fig. [5T] show an irreg- 
ular dependence on the vibrational state of the mol ecule . 
This has already been seen for H -f H2 coUisions [2l[lT3- 
In contrast, experimental measurements of quench ing 
rate coefficients for Cs -I- Cs2 collisions |l06l , llOTl j do 
not show any dependence on the vibrational state of the 
molecule. The differences between these systems can be 
explained based on the following considerations. 

First, the theoretical study applies to spin-polarized 
atom - dimer alkali-metal systems whereas it is not the 
case for the experiment. A full theoretical treatment 
should involve the electronic and nuclear spins of the al- 
kali metal atoms as well as couplings between electronic 
surfaces of different spins. This is beyond the scope of 
quantum dynamics calculations at present and will in- 
volve significant new code development and massive com- 
putational efforts. 

Second, the dynamics of the two systems are differ- 
ent. The lithium system is lighter and has a mo re at - 
tractive three-body term than the cesium system jl05| [. 
The triatomic adiabatic potential energy curves are well 
separated for a light system such as Lis and very dense 
for a heavy system such as CS3. The density of states 
strongly infiuences the nature of vibrational relaxation. 
The effect of the density of states has been illustrated in 
calculations for the Li -I- Li2 system which exclude the 
three-body term. Removing the three-body term makes 
the energy levels more sparse. This results in a more reg- 
ular and monotonic dependence of the rate coefficients on 
the vibrational levels ?; = 3 — 9 as illustrated in Figure [H] 
This conclus ion i s in agreement with a previous work of 
Bodo et al. |ll7{ for the H -f II2 system. When they 
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increased the density of states of the triatomic system, 
they found no significant dependence on the vibrational 
states of the molecule. However, for v = 10, they ob- 
tained the same results with and without the three-body 
term. The three-body term thus appears to be less signif- 
icant for collisions of molecules in high vibrational levels. 
Since the three-body term vanishes at large separations 
it is only important in the short-range interaction region 
while high vibrational states involve spatially extended 
molecules and sample the long range part of the interac- 
tion potential. Therefore, three-body terms, which are 
the most difficult interaction energy terms to compute 
numerically for triatomic systems, may be neglected to 
a first approximation in dynamics of highly vibrationally 
excited molecules. 

In the experiments, rate coefficients have been mea- 
sured fo r tem peratures of 40 x 10^^ K [l07l| and 60 x 
10~^ K jl06l | for Cs + Cs2 while the theoretical rate co- 
efficients for Li + Li2 were reported for a temperature 
of 10~^ K which corresponds to the Wigner threshold 
regime. Thus, it is likely that the experimental mea- 
surements did not probe the limiting values of the rate 
coefficients in the Wigner regime. 



4- Reactions of heteronuclear and isotopically substituted 
alkali-metal dimer systems 

Reactive collisions involving heteronuclear molecules 
are currently of great interest. A major goal of recent ex- 
periments has been to produce heteronuclear alkali-metal 
dimers in their e l ectro nic gro und s t ate. Examp les includ e 
RbCs [ll3, [M [Hi, NaCs [HO, [Hp, K Rb fl2l [HI, 
LiCs (124| . and mixed isotopes ^Li'^Li jl25l |. ^Quan- 
tum dynamics calculations for heteronuclear systems 
are more difficult. They are more interesting from a 
chemistry perspective because the reactive channels in 
collisions of heteronuclear molecules can be distinguished 
from inelastic channels. 

Cvitas et al. [H, [93] have explored the quantum dy- 
namics of ^Li -I- ^LfU{v = 0,j = 0), '^Li -I- ^U2{v = 
0,j = 1), + 7Li2(u = 0, j = 0), and ^Li + ^Li^U{v = 
0,j = 0) collisions. The J — cross sections for elas- 
tic scattering and chemical reactions in collisions of ^Li 
with ^Li^Li(w — 0,j = 0) are presented in the upper 
panel in Fig. [22] for a wide range of collision energies. 
The cross sections are believed to be converged for en- 
ergies up to 10""* K. The reactive process which leads 
to ^Li -I- ^Li2(u — 0, j = 0) dominates over the elastic 
scattering at ultralow energies. However, the reactive 
process is less efficient than the vibrational relaxation 
process in homonuclear alkali systems discussed above. 
For instance, at a collision energy of 10~^ K, the ratio 
of the cross sections for reactive and elastic collisions is 
about two orders of magnitude for the heteronuclear sys- 
tem presented in Fig. I22l whereas it is about three orders 
of magnitude for homonuclear systems. Cvitas et al. at- 
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FIG. 22: Elastic and reactive s-wave cross sections for ^Li 4- 
^CUiv = 0,j = 0) (upper panel) and '^Li -I- %i2{v = 0, j = 
1) (lower panel). Reproduced with permission from Cvitas et 
al. [ill. 



tributed the smaller ratio to the presence of only one 
open channel in the heteronuclear reaction. 

The J — 1 elastic and reactive cross sections for ^Li 
-I- ^Li2(u — Q,j — 1) collisions are presented in the lower 
panel of Fig. [521 The reactive process leading to the 
formation of ^Li^Li molecules is slightly more efficient 
than elastic scattering in the Wigner regime. The other 
possible collision processes are ^Li -I- ^Li^Li(?; = 0,j = 0) 
and ^Li + ^Li2(f = 0,j = 0). However, only elastic 
scattering occurs in these systems at ultralow collisions 
energies. 

These results provide important implications for the 
experiments on the production of ^Li''Li in an ultracold 
mixture of ^Li and ^Li atoms. Cvitas et al. proposed to 
remove quickly the ^Li atomic gas after the formation 
of the ^Li^Li molecules in their ground state to prevent 
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the destructive reactive processes. But keeping the ^Li 
atomic gas is recommended for sympathetic coohng 
of the ^Li^Li molecules since only elastic collisions 
are possible. Removing the ^Li atoms will leave the 
fermionic heteronuclear ^Li^'Li molecules in the trap. 
In the absence of ^Li atoms, evaporative cooling by 
collisions between fermionic ^Li^'Li molecules will not 
be effective since s-wave collisions of identical fermionic 
dimers will be suppressed due to the Pauli exclusion 
principle. 

Table IIVI provides a compilation of zero-temperature 
quenching rate coefficients for different alkali-metal 
trimer systems. 



IV. INELASTIC MOLECULE - MOLECULE 
COLLISIONS 

Many of the studies of cold and ultracold molecules 
have focused on reactive and non-reactive scattering in 
atom - molecule collisions. At high densities of trapped 
molecules, molecule - molecule collisions need to be con- 
sidered. The presence of rotational and vibrational de- 
grees of freedom in both collision partners make molecule 
- molecule systems especially interesting. However, quan- 
tum dynamics calculations of molecule - molecule colli- 
sions are significantly more challenging. Most dynamical 
calculations of molecule - molecule scattering have relied 
on the rigid rotor approximation and some recent studies 
have adopted the coupled-states approximation. Calcula- 
tions performed at high collision energies have employed 
more approximate methods based on semi-classical tech- 
niques. Here, we present a brief account of recent dy- 
namical calculations for the H2 + H2 system and discuss 
some ongoing work on full-dimensional quantum calcu- 
lations of ro- vibrational transitions in H2— H2 collisions. 
A brief discussion of hyperfine transitions in molecule - 
molecule systems is also provided for the illustrative ex- 
amples of the O2 + O2 and OH + OH / OD + OD 
systems. 



A. Molecules in the ground vibrational state 

The H2 + H2 system is the simplest neutral tetra- 
atomic system and it serves as a prototype for describing 
collisions between diatomic molecules. Though there 
have been a number of experimental and theoretical 
studies over the pas t several years on the H2 + H2 
system (see Ref. |l26l | and references therein), only a few 
studies have explored the collision dynamics in the cold 
and ultracold regime. 

Forrey |l27| presented a study of rotational tran- 
sitions in H2(w = 0,j = 2) + H2(u' = 0,j' = 2) 
collisions at cold and ultracold collision energies using 
a rigid rotor model. He calculated the real and imagi- 
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TABLE IV: Zero-temperature quenching rate coefficients for 
different atom - dimer alkali metal systems. 



nary components of the complex scattering length for 
H2(u = 0,j = 2,4,6,8) -I- H2(w' = 0,j' = j) collisions 
and showed that the imaginary parts decrease with 
increasing The imaginary parts of the scattering 

lengths are found to be small compared to the real 
parts, leading to small inelastic cross sections. Mate 
et al. (l28i | reported an experimental study of the rate 
coefficient for the H2(w = 0, j = 0) H2(w' = 0,j' = 0) 
}i2{vf = 0,jf = 0) -I- H2(t)^ = 0,j'f = 2) transition 
at temperatures between 2 K and 110 K. They found 
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FIG. 23: Elastic cross section of n2{v = 0, j = 0) + R2{v' = 
0; j' = 0) as a function of the collision energy. Reproduced 
with permission from Lee et al. ^126] . 



good agreement between the experimental results and 
quantum dynamics calculations based on the rigid rotor 
mode l usi ng a PES repor ted by Diep and Johnson 
(DJ) [l29l |. Montero et al. jl30| [ have investigated cold 
inelastic collisions of n-H2 molecules in the ground 
vibrational state, using a 3:1 gas mixture of ortho and 
para hydrogen. They obtained good agreement between 
the experimental data and the theoretical calculations 
based on the DJ PES. 

Using a quant um f ormalism based on the rigid rotor 
model, Lee et al. jl26| have recently presented a compar- 
ative analysis of cross sections for rotationally-inelastic 
collisions between H2 molecules at low and ultra-low 
energies. The elastic cross sections for H2(u = 0,j — 0) 
-I- H2(w' = 0,j' = 0) collisions obtained in this study 
are presented in Fig. [23] for two different PESs. The 
limiting value of the elastic scattering cross section in 
the ultra-low energy regime is 1.91 x 10"^'^ cm^ with 
the PES of B oothroyd, Martin, Keogh and Peterson 
(BMKP) [ml and 1.74 x 10"" cm^ with the DJ 
PES [l29| . For low collision energies the dynamics 
is sensitive to higher-order anisotropic terms in the 
angular expansion of the interaction potential. A diatom 
- diatom scattering length of 5.88 A was obtained for 
the DJ PES and 6.16 A for the BMKP PES. Cross 
sections for the quenching of rotationally excited H2 
molecules in H2(w = 0, j = 2) + H2(w' — 0,j' — 0) and 
H2(?; = 0, j = 2) -I- H2(w' = 0,/ = 2) collisions were 
presented at low and ultra-low energies. 

Quantum calculations of rotational relaxation of CO 
in cold and ultracold collisions with H2 have recently 
been performed by Yang et al. |l32l Il33j |. They 
reported quenching rate coefhcients for j = 1 — 3 of 
the CO molecule in collisions with both ortho- and 
para-H2 



cross sections exhibit a number of narrow resonances 
for collision energies between 1.0 — 40.0 cm~^. The 
signatures of these resonances are present in the tem- 
perature dependence of the rate coefhcient which shows 
broad oscillat ory features in the temperature range of 
10^2 „ 50 K p^ . 

Bohn and co-workers have reported extensive calcu- 
lations of hyperfine transitions in ultracold molecule - 
molecule collisions. Avdeenkov and Bohn studied ultra- 
cold collisions between O2 molecules 13^. They used a 
quantum mechanical formalism based on the rigid rotor 
model and included the electronic spin structure of the 
O2 molecules. Couplings between the rotational angu- 
lar momentum and the electronic spin of the molecules 
lead to rotational fine structure. Avdeenkov and Bohn 
discussed the elastic and inelastic loss processes in ^''02 
+ "O2 and 1^02 + ^^02 colhsions. They found that 
for collision energies below 10"^ K, elastic collision cross 
sections are larger than inelastic spin-flipping transi- 
tions. Based on relative magnitudes of elastic and in- 
elastic spin-flipping cross sections they concluded that 
^^02 molecules would be good candidates for evaporative 
cooling. In contrast, for ^^02 -I- ^^02 collisions, inelastic 
processes are more efficient than elastic collisions so that 
^^02 molecules are prone to collisional trap loss. 




132|. Due to the relatively deep van der 



Waals interaction potential for the H2— CO system the 



Energy(K) 



FIG. 24: Elastic and inelastic cross sections for OD + OD 
(thick gray line) and OH + OH (thin black line) collisions for 
an applied electric field of e=100 V/cm. Solid and dashed 
lines refer to elastic and inelastic cross sections, res pectiv ely. 
Reproduced with permission from Avdeenkov et al. [l37j . 



Avdeenkov and Boh n also studied ul traco ld collisions 
between OH [Ullilil and OD radicals in the pres- 
ence of an applied electric field. They showed that jl37| , 
elastic scattering is more efficient than inelastic processes 
for ultracold collisions between fermionic OD molecules, 
inhibiting state changing collisions. The energy depen- 
dence of elastic and inelastic cross sections for OH + OH 
and OD -I- OD collisions is illustrated in Fig. [24] for an 
applied electric field of e=100 V/cm. While the elastic 
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cross sections approach finite values in the Wigner regime 
for the bosonic system of OH molecules and the fermionic 
system of OD molecules, the inelastic cross sections ex- 
hibit a totally different behavior. At ultra-low energies 
when molecules interact with an electric field, s-wave 
scattering yields an inelastic cross section diverging as 

— 1/2 

^coii ^'^^ bosonic systems, while p-wave scattering yields 
an inelastic cross section vanishing as E^JJ^^ for fermionic 
systems. Thus, the inelastic processes are suppressed for 
the fermionic system. The differences are attributed to 
the bosonic/fermionic character of the molecules and to 
the applied electric field. In the absence of an applied 
electric field, the elastic cross section of the fermionic sys- 
tem would decrease as E^^;; , faster than the inelastic cross 

section, which decreases as E^^j^j . Ticknor and Bohn (l38| 
subsequently studied OH— OH collisions in the presence 
of a magnetic field. They showed that magnetic fields 
of several thousand gauss reduce inelastic collisions by 
about two orders of magnitude. Based on these results, 
they concluded that magnetic trapping may be favorable 
for OH molecules. 



B. Vibrationally inelastic transitions 

The theoretical studies presented above for molecule - 
molecule collisi ons refer to rigid rotor molecules. Pogreb- 
nya and Clary |l39l . Il40l | have investigated vibrational re- 
laxation in collisions of hydrogen molecules using a full- 
dimensional quantum dynamics formalism with an an- 
gular momentum decoupling approximation in the body 
fixed frame. They studied H2(w = l,j) -|- ii2{v = 0,j') 
collisions involving both para-para and ortho-ortho com- 
binations for collision energies of 1 meV (11.6 K) to 1 eV 
(11604 K). Time-dependent quantum mechanical calcu- 
lations based on the multiconfiguration time-dependent 
Hartree approach have been recently applied to study 
full-d imensiona l quantum dynamics of the H2— H2 sys- 
tem [l4ll . Il42l |. These methods are suitable at high 
collision energies and have not been applied to cold 
and ultracold collisions. Very recently, vibrational en- 
ergy transfer mechanism for ultracold collisions between 
para- h ydro gen molecules has been explored by Quemener 
et al. |l43l l using a full-dimensional quantum theoret- 
ical for malis m implemented in a new code written by 
Krems [l44j |. The quantum scattering calculations are 
based on the theory d escri bed by Arthurs and Dal- 
garno Takayanagi [l46| . Green |l47l |. and Alexan- 

der and DePristo [148i] . The H2 molecules were initially 
in different quantum states characterized by the vibra- 
tional quantum number v and the rotational angular mo- 
mentum j. A combination of two ro- vibrational states of 
H2 was referred to as a combined molecular state (CMS). 
A CMS, denoted as {vjv'j'), represents a unique quan- 
tum state of the diatom - diatom system before or after 
a collision. 

The cross sections for H2(w = l,j = 0) -I- H2(w — 



H,(v=1,j) + H,(v=0,j') 




Final combined molecular state 



FIG. 25; Elastic and inelastic cross sections for the collisions 
H2(« = l,j = 0) + H2(^ = 0,j = 0), H2(« = l,j = 2) + 
H2(« = 0, j = 0), and Hziv = 1, j = 0) + Hziv = 0, j = 2): 
total cross sections (upper panel); state-to-state cross sec- 
tions at 10~^ K (lower panel). Adapted with permission from 
Quemener et al. ,143.] . 



0,j = 0), H2(w = l,j = 2) + H2(f - 0,j = 0) and 
H2(w = 1, j = 0) + ll2{v = 0, j — 2) colhsions are pre- 
sented in the upper panel of Fig. [55] Using the CMS 
notation, this corresponds, respectively, to initial CMSs 
(1000), (1200) and (1002). The elastic cross sections are 
almost independent of the different initial ro-vibrational 
states of the molecules, but the inelastic cross sections are 
strongly dependent on the initial rotational and vibra- 
tional levels of the H2 molecules. At 10 ~^ K, the inelastic 
relaxation of H2 (u = 1, j — 0) is almost six orders of mag- 
nitude more efficient in collisions of H2(w = 0, j — 2) than 
in collisions with H2(u = 0, j = 0) and two orders of mag- 
nitude more efficient than in collisions of H2(w = 1, j = 2) 
with H2(w = 0,j = 0). At 25.45 K, which corresponds 
to the energy difference between the CMSs (1002) and 
(1200), the inelastic cross sections become comparable for 
U^{v = l,j = 2) -f H2(z; = 0, J = 0) and li2{v = 1, j = 0) 
+ H2(u — 0,j — 2) collisions. The inelastic scatter- 
ing depends on the type and the combination of ro- 
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vibrational levels involved in the collision, whether it in- 
volves ground state molecules, H2(w = 0,j = 0), vibra- 
tionally excited molecules, H2(w = 1, j = 0), rotationally 
excited molecules H2(w = 0,j = 2), or ro-vibrationally 
excited molecules, H2(w — l,j — 2). H2 molecules 
are weakly interacting and characterized by shallow van 
der Waals interaction at large separations and the com- 
puted rate coefBcients are not representative of strongly 
interacting alkali- metal dimer systems. For example, 
Mukaiyama et al. [l09j | reported an inelastic rate coef- 
ficient of 5.1 X 10~^^ cm^ s"'^ for collisions between two 
weakly bound Na2 molecules created by the Fesh bach 
resonance method. In a similar study Syassen et al. |llClj | 
reported a rate coefficient of 3 x 10^^" cm'^ s^^ for col- 
lisions between two Rb2 molecules. Zahzam et al. jl07j | 
estimated a rate coefficient of about 10~^^ cm^ s~ ^ for 
collisions between Cs2 molecules. Ferlaino et al. (l49j | 
measured a rate coefficient of 9 x 10~^^ cm^ for colli- 
sions between Cs2 molecules in the non-halo regime. The 
latter authors also measured the variation of the inelastic 
rate coefficient as a function of the atom - atom scatter- 
ing length for collisions between tunable halo dimers of 
Cs2. The large values of the inelastic rate coefficients 
for alkali-metal dimer systems are attributed to strong 
inelastic couplings and deeper potential energy wells. 

The state-to-state cross sections for three ro- 
vibrational combinations of the H2— H2 system are 
presented in the lower panel of Fig. 1251 The magnitude 
of the inelastic cross sections depends on the propensity 
of the diatom - diatom system to conserve the internal 
energy and the total ro tatio nal angular momentum of 
the colliding molecules jl43l |. The final state-to-statc 
distribution in H2(w = IJ = 0) + li2{v = 0,j = 0) 
collisions shows that there is no preferential popula- 
tion of rotational quantum number of either of the 
colliding molecules. The conservation of the total 
rotational angular momentum would entail a large 
change in the internal energy of the molecules. Thus, 
the purely vibrational transition (1000) (0000) is 
not efficient because the energy gap is large. On the 
other hand, (near) conservation of the internal energy 
requires a large change in the total rotational angular 
momentum of the colliding molecules and the transition 
(1000) (0800) is not dominant either. However, the 
state-to-state cross sections for H2(w — l,j — 2) -|- 
112 {v = 0,j = 0) collisions indicate that the transition 
(1200) — >■ (1000) is more efficient than all the other 
transitions combined. For this transition, the total 
rotational angular momentum change and the internal 
energy transfer are both minimized, leading to a more 
efficient energy transfer process. The state-to-state 
cross sections for H2(w = 1, j = 0) -I- H2(w = 0,j = 2) 
collisions present an interesting scenario in which the 
transition (1002) (1200) is highly cflfcient and 
selective. In this case, the total rotational angular mo- 
mentum is conserved and the internal energy is almost 
unchanged (the energy gap between (1200) and (1002) 
is only 25.45 K). This creates a favorable situation 



leading to a near-resonant energy transfer process. 
This particular mechanism cannot occur in atom - 
diatom systems because simultaneous conservation of 
rotational angular momentum and internal energy of the 
molecule cannot occur. The mechanism is reminiscent of 
quasi-resonant energy transfer in collisions of rotation- 
ally excited diatomic molecules with atoms discussed 
previously i53i, 113, but with a purely quantum 
origin. The near-resonant process may be an important 
mechanism for collisional energy transfer in ultracold 
molecules formed by photoassociation of ultracold atoms 
and for chemical reactions producing identical molecules. 

Table |V] provides a compilation of zero-temperature 
rate coefficients for rotational and vibrational quenching 
in molecule - molecule systems. 
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TABLE V: Zero-temperature inelastic rate coefficients for dif- 
ferent molecule - molecule systems. 
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V. SUMMARY AND OUTLOOK 

In this chapter we have given an overview of recent 
theoretical studies of atom - molecule and molecule - 
molecule collisions at cold and ultracold temperatures. 
Though such systems have been extensively studied 
at higher collision energies over the last few decades, 
the new experimental breakthroughs in creating dense 
samples of cold and ultracold molecules have provided 
unprecedented opportunities to explore elastic, inelastic, 
and reactive collisions at temperatures close to absolute 
zero. These studies have revealed unique aspects of 
molecular collisions and energy transfer mechanisms that 
are otherwise not evident in thermal energy collisions. 

The long duration of collisions combined with large de 
Broglie wavelengths at cold and ultracold temperatures 
leads to interesting quantum effects. Calculations 
have shown that reactions with insurmountable energy 
barriers may still occur at temperatures close to absolute 
zero, and in certain cases, with appreciable rate coeffi- 
cients. Such tunneling dominated reactions have been 
the topic of many recent investigations and may soon be 
amenable to experimental investigation in the cold and 
ultracold regime. That the rates of these reactions can 
be enhanced by vibrational excitation of the molecule 
is an interesting scenario for experimental studies of 
ultracold chemical reactions. 

The studies of alkali-metal homonuclear and het- 
eronuclear trimer systems have stimulated considerable 
experimental interest in investigating chemical reactivity 
at ultracold temperatures. The challenge for theory 
is to describe collisions involving highly vibrationally 
excited molecules. Recent theoretical calculations have 
indicated that the three-body interaction potential can 
be neglected for highly vibrationally excited molecules 



offering significant savings in computational effort. Even 
so, heavier alkali-metal trimer systems pose a daunting 
computational challenge. 

The quantitative description of ultracold molecule - 
molecule collisions is another challenging topic. The 
recent progress on the H2— H2 system will be difficult 
to implement for heavier systems due to the large 
number of ro-vibrational levels of the molecules. The 
study of H2— H2 collisions has shown that, for certain 
combinations of ro-vibrational levels, the energy transfer 
may occur to specific final ro-vibrational states. In such 
cases, the calculations can use a much smaller basis set 
without compromising the accuracy. 

Currently there is substantial interest in controlling 
the collisional outcome using external electric and 
magnetic fields. While the idea of coherent control of 
molecular collisions and chemical reactivity has existed 
for a long time and some important progress has been 
achieved, the possibility of creating coherent and dense 
samples of molecules in specific quantum states has 
given further impetus to the field of controlled chemistry. 
We expect that the coming years will see a far greater 
activity in this direction driven by cold and ultracold 
molecules and also by the possibility of controlling 
chemical reactivity using external electric and magnetic 
fields. Electronically non-adiabatic effects in ultracold 
collisions is a largely unexplored area, which can be 
expected to attract much attention. 
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